On the locus formed by the maximum heights of projectile motion with air resistance 
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We present an analysis on the geometrical place formed by the set of maxima of the trajectories 
of a projectile launched in a media with linear drag. Such a place, the locus of apexes, is written 
in term of the Lambert W function in polar coordinates, confirming the special role played by this 
function in the problem. In order to characterize the locus, a study of its curvature is presented in 
two parameterizations, in terms of the launch angle and in the polar one. The angles of maximum 
curvature are compared with other important angles in the projectile problem. As an addendum, 
we find that the synchronous curve in this problem is a circle as in the drag-free case. 
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I. INTRODUCTION 

An amazing characteristic of some old fashion prob- 
lems is their endurement. The projectile motion is one 
of them. Being one of the main problems used to teach 
elementary physics, variations and not well known facts 
about it appear in the physics literature of the XXI cen- 
tury. A search in the web [l| or in the Science Citation 
Index gives an idea of this fact. Some of the recent studies 
deal with the problem of air resistance in the projectile 
motion and its pedagogical character made of it an ex- 
cellent example to introduce the Lambert W function, a 
special function. [2I] The Lambert W function is involved 
in many problems of interest for physicist and engineers, 
from the solution of the jet fuel problem to epidemicsfl] 
or, even. Helium atom eigenfunctions.0 One of those 
problems is the solution for the range R in the case that 
the air resistance has the from f = —mbv.^,^ 

In this paper we analyze the not well known fact of the 
geometrical place formed by the maxima of all the projec- 
tile trajectories at launch angle a and in the presence of 
a drag force proportional to the velocity, we shall denote 
this locus as Cm{s). The resulting locus becomes a Lam- 
bert W function of the polar coordinate r{9) departing 
from the origin. This problem raises as a natural contin- 
uation from the nice fact that in the drag-free case such 
a locus is an ellipse[^-Q with an universal eccentricity 
e = 73/2.1] 

The paper is organized as follows. In section [Hi the 
set of maxima for projectile trajectories moving under in 
the presence of air resistance is presented. In section Hill 
we find a closed form, in polar coordinates, to express 
such a geometrical place. Cm- In section ITVl we present 
a numerical calculation of the curvature of Cm using the 
polar angle and the launch angle as parameterizations. 
Additionally, we demonstrate that the synchronous curve 
is a circle as in the drag- free case in section IVl In section 
we conclude. 



II. THE PROJECTILE PROBLEM WITH AIR 
RESISTANCE 

Several approximations in order to consider the air re- 
sistance exist in the literature, the simplest is the linear 
case. In such a case the force is given by 



F = —mhv — mg], 



(1) 



where m is the mass of the projectile and b is the drag 
coefficient. The units of b are s~^. The velocity compo- 
nents are labeled as v = ui + wj, with u = dx/dt and 
w — dy/dt. The solutions for the position and velocity 
are obtained trough direct integration of Eq. ([T]) yielding 



x{t) 



if [1 



exp{~bt)] 



y(t) ^ !£o+£Z^[i_exp(-M)]-5t/6, 
for the coordinates, and 



u{t) — uoexp(— M), 
w{t) = (wo + g/b) exp(- 



bt) - g/b, 



(2) 



(3) 



(4) 
(5) 



for the speeds. We used the initial conditions x{0) = 
2/(0) = and uq = Vocosa and wq = Vbsina. Noticing 
that the terminal speed is g/b in the y axis. 

For the same initial speed Vq these solutions are func- 
tion of the launch angle a and the locus formed by the 
apexes is obtained if time is eliminated between the so- 
lutions in time in Eqs. ([2]) and ([3]), giving the equation 



Wq + g/b ^ g_ 
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and considering the value at the maximum, via dy/dx 
0. The corresponding solution is 
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esina — ln(l -t- esina), 



(8) 
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FIG. 1: (color online) Locus Cm(e) formed by the apexes of 
all the projectile trajectories (continuous line in blue) given 
by Eqs. ^ and ((8]l in rectangular coordinates or by Eq. (|13|l . 
the last one express Cm in polar coordinates and in term of 
the Lambert W function. The dashed red line is the ellipse 
of eccentricity e = VS/2 which represents the drag-free case, 
i.e. Cm(0). The parameters are Vb = 10 and e = 0.1. 



where we introduce the dimensionless perturbative pa- 
rameter £ = bVo/g, the dimensionless length p = Vq /g, 
and noticing that y can be expressed as ^. An alter- 
native procedure consists in set the derivative dy/dt to 
zero to obtain the time of flight to the apex of the tra- 
jectory and, evaluate the coordinates at that time. The 
points {xm,ym) conform the locus of apexes Cmi^) for 
all parabolic trajectories as a function of the launch an- 
gle a. In Fig. [1] we plot Cm{e) described by Eqs. ([7]) 
and ([5]), for the drag-free case (in dashed red line) and 
for £ = 0.1 in continuous blue line. Several projectile 
trajectories are plotted in thin black lines. The locus of 
apexes Cm{£) defined by {xm,ym) of Eqs. ([7]) and ^ 
is described parametrically by the launch angle a and it 
changes for different values of £ = bVo/g. In the next 
section we shall find a description of Cmis) in terms of 
polar coordinates and in a closed form using the Lambert 
W function. 



rearranging terms it must be expressed as 



cos 61™ exp (^-£2 sin 
— cosasinQ;exp(— esina). 



(9) 



The Ihs depends on r„i and 9m meanwhile the rhs de- 
pends on a, however the last angle is a function of 9^ 
and reads as 



tan 9 J, 



1 (£sinQ; — ln(l + £sinQ;)) 



cos a sm a 



(10) 



l-\-s sin a 



by making tan0,„ — ym/^m from Eqs. ([7]) and 
In order to obtain r{9„i) = r„ii9m)/p we set 

f{a{9m)) = cosQ!sinQ;cxp(— £sinQ;), (11) 

since Eq. pH)) allows us to have, implicitly, a{9m)- We 
shall return to this point later. Hence, we can write Eq. 
& as 



— £ sm ( 



r{9m) exp (-£ sm9„i f{9,n)) = - 



= — £ tan( 



m /(a). 
(12) 

Where we multiplied both sides of Eq. ^ by 
— £^sin6'm. Setting z = — £^tan6',„ f{a) and W{z) = 
—e^sm9m r{9m) in Eq. (fT2| . it shall have the familiar 
Lambert W function form, z = W{z) ex-p{W{z)), from 
which we can obtain f as 



1 



-W{- 



-e^ tan( 



/(«))■ 



(13) 



It is important to note that the argument of the Lam- 
bert function in this equation is negative for all the values 
£ > 0. W{x) remains real in the range x £ [— 1/e, 0) and 
have the branches denoted by and — 1.0 We select the 
principal branch, 0, since it is the bounded one, however, 
for values of £ > 1.1 there is a precision problem since the 
required argument values are near to —1/e = — exp(— 1). 
It is important to stress that in Eq. ([T^ the independent 
variable is the angle 9 and, it constitutes the parameter- 
ization of the curve Cm ■ 

The polar expression of Cm can also be written in terms 
of the tree function T{z) — —W{—z), giving 



f{9n 



-Tie^ i(,n9mf{a{9m)))- 



(14) 



III. THE LOCUS Cm AS A LAMBERT W 
FUNCTION 



In order to obtain an analytical closed form of the lo- 
cus we change the variables to polar ones, i.e., Xm — 
COS 9m and ym — sin 9m • 

The selection of a de- 
scription departing from that origin instead of the center 
or the focus of the ellipse is because the resulting geomet- 
rical place is no longer symmetric and the only invariant 
point is just the launching origin. We substitute the po- 
lar forms of Xm and ym into equations ([7]) and ([5]) and 



We recover the drag-free result 

sin 9ri 



' 1 + 3 sin^ 9r, 



(15) 



when £ — !■ 0. An explanation of this unfamiliar form of 
an ellipse is given in appendix El followed by a discussion 
about the £ — ^ limit of expression ([T^ in appendix |B] 
Formula ([T^ exhibits the deep relationship between 
the Lambert W function and the linear drag force pro- 
jectile problem, since not only the range is given as this 
function HQ. The problem open the opportunity to 
study the W function in polar coordinates, that, almost 
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FIG. 2: (color online) The angle S as a function of the launch 
angle a for two different values of e and their inverses. 



in the review of reference! [2:|, is absent. Even when it is 
possible to write the locus in terms of ym{xm) this form 
does not shows the formal elegance of relation ([T^ . 

Now we return to equation (jlOp since we need to solve 
explicitly it in order to have the function a{9m)- This 
task is not trivial since even when we approximate the 
rhs in expression (1101) up to first order in e, 



tana e . 
tanWrn ~ — - — (1 + - sma) 



(16) 



the inversion is not easy. A way to do the inversion is 
to expand in a Taylor series the rhs and then invert the 
series term by term.Q Using Mathematica to perform 
this procedure up to 0(18), we obtain as a result 



a{9) « arctan(2tan6') - ie(2tan6')2 
+-tj(27£ - 10£3)(2tan6i))4 + 



'(2tan( 



(17) 



The e-independent terms had been resumated to yield 
arctan(2tan6'). However, the series does not converge 
for values in the argument larger than 1. The reason is 
the small convergence ratio for the Taylor expansion of 
arcsin(-). 

An easier way to perform the inversion is to evaluate 
9jn (a) using Eq. ([TU)) and plot the points (6*™ (a) , a) , the 
result is shown in figure [51 The result is in agreement 
with the plot of Eq. p7)) up to its convergence ratio 
and it is not shown. Notice that this method is exact in 
the sense that we can obtain as many pair of numbers 
as we need, a function is, finally, a relation one to one 
between two sets of real numbers. Another result is to 
obtain the derivative da/d9, since it shall be needed in 
the following sections. To this end, we note that both 
functions increase monotonically and their derivatives are 
not zero, except at the interval end. Hence, we can use 
the inverse function theorem in order to obtain 



da 
'd9 



1 



da 



(18) 



FIG. 3: (color online) (a) First and (b) second derivatives of 
a as function of 9 for various values of parameter e. Note that 
major changes occur for 9 < 7r/4. 



The result is shown in Fig. Ela) as well as the second 
derivative in Fig. EJb). The second derivative is calcu- 
lated using an approximation to the slope to the func- 
tion previously calculated and using 10000 points in the 
interval [0, 7r/2]. A smaller number of points could be 
considered. 



IV. THE CURVATURE OF C^. 
A. Polar angle parameterization. 

In the drag-free situation. Cm is an ellipse and its de- 
scription is well know, however, in the presence of linear 
drag this is not the case. We do not expect that the lo- 
cus could be a conic section and henceforth we need to 
characterize it. It is usual to consider curvature, radius 
of curvature or the length of arc in order to characterize 
a locus. In the present case we consider the curvature 
of Cm in both parameterizations, first with the polar an- 
gle 9m and secondly with the launch angle a. We left the 
calculus of the length of arc to a posterior work, since the 
calculations became increasingly complex and the goal of 
the present section is to start the understanding of Cm 
and to illustrate the way it can be done using the Lam- 
bert W function. Here and in the rest of the section we 
drop, for clearness, the subindex m in fm and 9m- 

The corresponding formula for the curvature K for po- 
lar coordinates is Oil 



K 



1 + 2f 



?2 



rroo 



P (f2 -I- f2)3/2 



(19) 



in order to use Eq. ([T3|) . Here the subindex 9 corresponds 
a derivative respect to that variable. 

A direct calculation on the drag- free f{9) of Eq. p^ 



yields to 



1 



(5 - 3cos26i)S 



2^/2p (1 + 3 sin^ 61)3 (47 - 60 cos 26* + 21 cos 461)3/2 ' 

(20) 

which have a maximum at 6* = 1/2 arctan(4/3) « 0.464. 
A graph of this results appears in figure red line). Note 
that this value is different from that we obtain if we eval- 
uate 9(a — 7r/4) = 1/2, the launch angle of maximum 
range. The maximum curvature happens at a smaller 
angle than the angle of maximun range. It is interesting 
to note that the angle 26* = arctan(4/3) corresponds to 
a triangle which sides fulfill the relation 3^ + 4^ = 5^, a 
Pythagoras' triple. 

Using the numerical results for a{9) from the previous 
section, it is possible to carry on the calculation of K 
(see figure |4]) performing the derivatives of f from Eq. 
(|13p in a direct form and evaluating numerically the re- 
quired values of a and its derivatives. For the required 
derivatives oi W{z) we used the expressions 



W{x) 



and 



d^W{x) - exp(-2Vl^(a;))(W^(a;) + 2) 



(1 + T4^(a;))3 



(21) 



(22) 



Using this method, we obtained good results for values of 
£ up to ~ 1 but we require to calculate arguments of the 
Lambert W function near the limit z — — 1/e for larger 
values of e. The reliability of our numerical result was 
done comparing the first and second derivatives of r(9) 
with those corresponding to the ellipse. 

As can be seen in figure |4l K present a maximum in 
all the cases which can be calculated as well. We left to 
an ulterior work the analysis of the maxima distribution 
as a function of the pertubative parameter, that is not 
the case for the curvature with a parameterization as we 
shall see in the next section. 



B. Launch angle parameterization. 

For the launch angle parameterization of Cm we shall 
use expression [l3| 



K = 



x'y" - y'x" 
(x'2 + y'2)3/2 



(23) 



for calculate the curvature from the rectangular form of 
equations (O and ([8]) , where the primes denote derivative 
respect the parameterization variable, a in this case. 
A direct calculation yields 
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FIG. 4: (color online) Curvature of Cm(e) using the polar 
angle as the parameter from Eq. (|19|) . In red line appears the 
corresponding result for the ellipse, Eq. (|20|) . The maximum 
happens At 6 = 1/2 arctan(4/3), different from the value of 
maximum range q = 7r/4. 



With 



and 



V2{a) 



16 + 6e2 - cos 2a + 2e2 cos 4a 
-|-30e sin a ~ e sin 3a + e sin 5a, 



5 + 3 cos 4a + Se^ - Ae'^ cos 2a+ 
£2 COS 4a + 10£ sin a — 5£ sin 3a- 
£ sin 5a. 



(25) 



(26) 



In the limit e 



Kq 



0, we recover the drag- free curvature 
16^/2 



1 



p (5 + 3cos4a)3/2' 



(27) 



which have a maximum at a = 7r/4 in the interval a £ 
[0,7r/2], as expected. A plot of pK{a) for several values 
of £ beginning at zero and ending at £ = 10 appears in 
Figureinja). In red appears the drag-free case. Note that 
both extremal values increase for increasing e value as 
Pk(0) w (16-h6£-6£2) and pK{n/2) sa e. Notice that for 
small £, the k crosses the drag- free curvature kq. 

The angles a* at which n attain their maxima are ob- 
tained in the usual way and requires to solve, numerical 
or graphically, the equation 



3(1 -H£ sin a*) sin 2a* x Qi{a*)Q2{a*)- 
£COsa*S3(a*)Q4(a*) =0, 



with 



Qi = 4 (3-}-£2)cos2a* +£(-4£- 15sina* 
5 sin 3a*); 



(28) 



(29) 



V2 Vlja) 



X (1 + £sina) 



(24) 



16 + 6£2 - 8£2 COS 2a* + 2e^ cos 4a* 
30£ sin a* — £ sin 3a* -I- £ sin 5a* ; 



(30) 
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FIG. 5: (color online) (a) Curvature as a function of the 
launch angle a for increasing values of parameter e according 
to equation (|24|l . The values of e are indicated in the inset, 
(b) Several important angles as a function of dimensionless 
parameter e are plotted. In lines and diamonds appears the 
angle, a* , at which the curvature is maximum. In red circles 
the angle at which the range is maximum according to exact 
solution, Eq. (1341) . and in blue crosses the same angle accord- 
ing to Eq. (1351) (see text for discussion) . In a dashed blue line 
the angle at which skewness is maximum is plotted. 



5 + 3 cos Aa* + Se^ - 46^ cos 2a* 

cos Aa* + lOe sin a* — 
5e sin 3a* + e sin 5a* ; 



(31) 



and 

24 



70 H 
2(5 



36e2 - 16(1 + 3e2)cos2a*- 
f 6 e^) cos 4a* + 154e sin a* 



3l£sin3 at32) 



"7e sin 5a* 



In Figure El^b) the calculated values of a* as a function 
of e appear. This angle is between the optimal angle for 
maximum range (red circles and blue crosses) and the 
angle for the greatest forward skew (dashed line) . 



^ske 



^((i?+/2)i/3 + (I?_/2)V3_2) 



where D± = ±3e\/3(27£2 - 32) + 27e^ 



(33) 

16, valid for 

e > a/32/27. In Figure [5][b) the optimal angles are 
drawn, in red circles the exact result in terms of Lambert 
W functionji, ^ 



amax,s = arcsm 



exp(W^(£^) + 1) - 1 
and the approximated result [3| 

W{eye) 



^max ,W 



(34) 



(35) 



Both expressions are equivalent for large e but differ at 
small e, as expected. 
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FIG. 6: (color online) Orbits launched at different angles and 
the corresponding geometrical place Cm(e) for e = 80. See 
text for explanation. 



Meanwhile the difference between these angles at small 
pertubative parameter is unimportant, at large £ the be- 
havior of the corresponding trajectories is different. One 
reason is the large asymmetry in the locus formed by the 
set of apexes. In Fig. [SJa) we plotted Cm{£) and the 
corresponding trajectories for the different launch angles 
for £ = 80. The blue line corresponds to Cmis), note that 
the maximum height is y ~ 0.12 in contrast to y ~ 5.0 
for the drag-free case, however, this can be the case of 
small friction parameter b and large initial velocity Vq 
giving a large £ value. In a black line appears the orbit 
launched at a*, in red line the corresponding to attain 
the maximum range and in blue dashed line the orbit 
with maximum skewness. 



V. THE SYNCHRONOUS CURVE 

In MacMillan's bookfzi] the calculation of the syn- 
chronous curve was done for the drag-free case. This 
curve is formed if many projectiles were fired simulta- 
neously from the same point, each one with at different 
launch angle and same initial speed Vq. The locus will be 
a circle of radius Vot and center in the point (0, —•^gt^), 
i.e. 



(36) 



Here, we shall demonstrate that a circle is the syn- 
chronous curve in the linear drag case as well. Following 
reference Q, we eliminate the launch angle a from the 
position solutions, in the present case they are equations 
and ([3]) . We write down cos a and sin a and rearrange 



6 



the terms to give, 
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COS a 



sm a 



Vo 1 - exp(-6) ' 
b y- -§z{l-bt^cxp{-bt)) 
Vo 1 - exp(-6i) 



(37) 



(38) 



Substituting these expressions in the identity cos^ a - 
sin^ a = 1 we obtain 



With 



= ^{l-bt-expi-bt)), 



the center and 



i?(t)--5^(l-exp(-6i)), 



(39) 



(40) 



(41) 



the radius. In order to recover the case where 6^0, 
we consider a Taylor expansion for the exponential up to 
second order in the exponential in Eq. (^0)) and up to 
first order in Eq. (|4T|) . The fact that this circle exists in 
the presence of a drag force is remarkable. 



VI. CONCLUSIONS 

We obtained an explicit form for the locus Cm com- 
posed by the set of maxima of all the trajectories of a 
projectile launched at an initial velocity Vq, and in the 
presence of a linear drag force, —mbv, i.e. Cm is the lo- 
cus of the apexes. In polar coordinates. Cm is written in 
terms of the principal branch of the Lambert W function 
for negative values. This represents the parameterization 
of the curve by the polar angle 9m only and gives Cm in 
a closed form and exhibits the deep relationship between 
the Lambert W function and the linear drag problem. 
The curvature of Cm was calculated for different values 
of the dimensionless parameter s = bVo /<? in two param- 
eterizations. The first one, the polar parameterization, 
shows a maximum that slightly departs from the drag- 
free case in 9 = 1/2 arctan(4/3). A wider exploration of 
the functional dependence respect to e is pending due to 
numerical accuracy in the calculation of the Lambert W 
function near the limit at a; = — 1/e. In the case of a pa- 
rameterization using the launch angle a there is not such 
a restriction. In this case, the curvature was calculated 
for a wide range of the parameter e yielding maximum 
at angle values larger than those corresponding to max- 
imum range. Comparison with the maximum skewness 
angle [l2| was also done and the difference is larger than 
the previous one. As an addendum, we demostrate that 
the synchronous curve, in this case, is a circle as in the 
drag-free case. 
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Appendix A: Polar form of an Ellipse with origin at 
bottom. 



Ellipse canonical form or the polar form with the origin 
considered in one of the focus are standard knowledge. 
In the present case, however, we require to consider the 
origin of coordinates located in the bottom of the ellipse, 
since, in the presence of a drag force, the launching origin 
is the only invariant point when we change the drag force 
value. To obtain the ellipse form, we depart from the 
drag-free solutions at the locus of the apexes. 



and 



Xm = psmacosa. 



P ■ 2 

ym = ^ sm a, 



(Al) 



(A2) 



being p = With the help of the trigonometric re- 
lations sin 2a = 2 sin a cos a and 2 sin^ = 1 — cos 2a we 
transform the upper equations into 



sin 2a 



cos 2a = 1 



P 



(A3) 
(A4) 



Taking the squares in both expressions, summing them 
and arranging terms, we arrive to 

— (—{l + 3smHm)~2sm9m]=0. (A5) 
P \ P J 

Where we used the polar coordinates Xm = fm cos 9m 
and j/„i — I'm sin 9m ■ The solutions are = and 



' my-'rn) ^ ^ ^ , 2 



(A6) 



1 + 3 sin"' 9m 

The second one is the required form for the ellipse. 

Appendix B: Drag- free limit for f 

In order to obtain the drag-free limit for the locus Cm 
given in Equation (|13p . wc note that 

(Bl) 



tan6' = (1/2) tana 



and that /(a) — sin a (^m) cos a (0m)- The first expres- 
sion is obtainable from the drag-free solutions Eqs. (|A1[) 
and (IA2p , and the second is obtained by setting 6 in 
Eq. (HI]). 
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The expansion in a power series of the Lambert W 
function up to first order is just the identity 0| and hence, 



hne. Using trigonometric identity sec^ a — tan^ a = 1 and 
Eq. (|Bip we obtain that 



= - 



1 

z'-^ sin 6jj 



sin a cos a. 



■T^(-e2 tan0„/(a)) 
■(-e^ tan e„rf{a)) 



(B2) 



cos^ 0„ 



= 1 + 3 sin^ 



(B3) 



= 2sin( 



Where we used relation (IBll) in order to obtain the last 



Using this result in the expression of r we obtain the 
desired result, Eq. ([T5|) . 
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